Auxiliary field method and analytical solutions of 
the Schrodinger equation with exponential 
potentials 

Bernard Silvestre-Brac 1 , Claude Semay 2 and Fabien 
Buisseret 2 

1 LPSC Univcrsitc Joseph Fourier, Grenoble 1, CNRS/IN2P3, Institut 
Polytcchnique de Grenoble, Avenue des Martyrs 53, F-38026 Grenoble-Cedex, 
France 

2 Groupe de Physique Nucleaire Theorique, Universite de Mons-Hainaut, 
Academie universitaire Wallonie-Bruxclles, Place du Pare 20, B-7000 Mons, 
Belgium 

E-mail: silvestre01psc.in2p3.fr, claude.semay9umh.ac.be, 
f abien.buisseretSumh. ac .be 

Abstract. The auxiliary field method is a new and efficient way to compute 
approximate analytical eigenenergies and eigenvectors of the Schrodinger 
equation. This method has already been successfully applied to the case of central 
potentials of power-law and logarithmic forms. In the present work, we show that 
the Schrodinger equation with exponential potentials of the form -«r* exp(— /3r) 
can also be analytically solved by using the auxiliary field method. Formulae 
giving the critical heights and the energy levels of these potentials are presented. 
Special attention is drawn on the Yukawa potential and the pure exponential one. 



PACS numbers: 03.65.Ge 



Submitted to: J. Phys. A: Math. Theor. 



Auxiliary field method for exponential potentials 



2 



1. Introduction 

Much work has been devoted to the computation of analytical solutions of the 
Schrodinger equation for the early years of quantum mechanics, especially when the 
Hamiltonian of the problem admits bound states. Indeed, it is often useful in physics 
to have a guess about the behavior of the eigenenergies of a system in terms of the 
various parameters of the Hamiltonian. This can be particularly relevant when one 
tries to fit the parameters of a model to some experimental data. 

There are many methods allowing to find approximate analytical solutions of the 
Schrodinger equation for bound state problems: WKB method, variational method, 
etc. (see pQ for example). In previous works [21 [3], we have proposed and extensively 
studied a new method to compute the eigenenergies and wave functions of a given 
Hamiltonian admitting bound states. This method is based on auxiliary fields (also 
known as einbein fields 14.) and can lead to analytical approximate results. Its general 
principle will be recalled in Section [2l Up to now, the auxiliary field method (AFM) 
has been mainly applied in the case of central potentials of the form ar x +br v . We have 
shown, in particular, that the qualitative features of the exact energy spectra are well 
reproduced and that the analytical formulae can be greatly improved by a comparison 
with a numerical resolution of the problem. Then they become very accurate [21 [3] . 

In the present work, we intend to apply the AFM to find approximate analytical 
formulae for central potentials of exponential form, that is — ar x e~@ r . We thus start 
with the following Hamiltonian 

H=^-ar x e- 0r , (1) 
2m 

where a and (3 are both positive and nonzero real parameters. It is interesting to 
define the new variables {x — (3r, q = p/(3} and to work with the dimensionless 
Hamiltonian H — 2mH/ 1 instead of H . Indeed, one has then 

H = q 2 -gx x e- x , (2) 



with 

2ma 



X+2 ' 



(3) 



The rewriting of the problem under this particular form emphasizes the fact that 
the dimensionless parameter g given by ^ contains all the relevant information of 
H. Only the eigenenergies e(<?) of ^ have thus to be computed; the "physical" 
eigenenergies will then be given by E(m,a,(3) — f3 2 e(g)/(2m), a relation which 
expresses the scaling properties of the Schrodinger equation. 

We first show in section [3] that the pure exponential potential (A = 0) can 
be analytically solved using the AFM provided that the energy spectrum of the 
Schrodinger equation with a linearly rising potential is analytically known. The pure 
exponential potential is common in hadronic physics for example, when a screened 
confining potential is considered [5J IB]- Then, we show that some general results can 
be obtained for A ^ in section[4l and write an explicit energy formula for the Yukawa 
potential (A = — I) in section[5l The importance of this potential in theoretical physics 
is such that it also deserves a specific study. We then present improved energy formulae 
for the pure exponential and the Yukawa potentials in section [BJ 

An interesting feature of exponential potentials is that they all admit a finite 
number of bound states that depends on the dimensionless parameter g, ruling the 



Auxiliary field method for exponential potentials 



3 



potential depth, and defined by (J3J. There exists thus "critical heights": Potential 
depths beyond which new bound states appear. We refer the reader to [7J for detailed 
explanations about how to compute critical heights in a given potential. Accurate 
formulae giving these critical heights are given for the pure exponential and the 
Yukawa potentials. Finally we sum up our results in section [7J The definition and the 
properties of the Lambert function, that will frequently appear in our calculations, 
are given in | Appendix A| Additional analytical approximations for the energy levels 
and some critical heights for the Yukawa potential are given in |Appcndix B| and 
|Appcndix C| respectively. 



2. Auxiliary field method 

2.1. General principle 

We recall here the main points of the method that we presented in [2J. The problem 
is to find an analytical approximate solution of the eigenequation H \*£) = E\ty), H 
being a Hamiltonian of standard radial form 

H = T(p 2 ) + V(r). (4) 

The procedure is the following. Let us assume that H A — T(p 2 ) + v P(r) admits 
bound states with an analytical spectrum, v being a real parameter. We thus know 
the solution of 

H A \y A (v))=E A {v)\y A {v)). (5) 
Let us define the function I(x) = K^ 1 {x) and the operator v from the relation 

The spectrum of the following Hamiltonian 

H{v) = T(p 2 ) + v P(r) + V -vP [I(y)\ , (7) 

is then analytically known provided that v - the auxiliary field - is assumed to be a 
real parameter. Using ([5]), we have indeed 

E{v)=E A (v)+V[I(v)]-vP[I{v)). (8) 

The physical spectrum is finally given by the eigenstates | v E'a(i / o)) and the 
eigenenergies E(vo), with vq minimizing the total energy ([8]), i.e. satisfying 

dE{v) 



dv 



= 0. (9) 

V = Vq 

This method can be justified as follows. The arbitrary potential V(r) is replaced 
by V(r, v) — v P(r) + V \I(y)\ — vP [1(f)] , a function of the auxiliary field v and of an 
analytically solvable potential P(r). The key point is that the auxiliary field should 
rigorously be eliminated as an operator from Hamiltonian ([7]). When doing this, one 

is led to 5 v V(r, v) = => = K(r), and to V(r, v) — V(r). This means that 

Hamiltonians (U|) and ([7J are equivalent up to an elimination of the auxiliary field as an 
operator, the auxiliary field being then defined by The AFM would consequently 
lead to the exact results if the real number v were replaced by the operator i*. As we 
use the real number v in order to have analytical solutions, the results are approximate. 
One has then u w (^^(i/q)] v \^ a (i/q)), with v and v given by ([6]) and (J9]) respectively 
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[2J. Our method can consequently be regarded as a "mean field approximation" with 
respect to a particular auxiliary field which is introduced to simplify the calculations. 

Let us notice that the method we have presented is not variational, and that it 
is exact when V(r) — P(r), with trivially v = 1. The main technical problem is the 
determination of analytical solutions for ((6|) and (J9]). Such a task has already been 
fulfilled in a satisfactory way for many interesting and nontrivial potentials [3] . As we 
will see in the latter, it can also be achieved for the Yukawa and the pure exponential 
potentials. 



2.2. Useful formulae 

We have shown in [2] that the eigenvalues of the Hamiltonian 

H ri (a) = ^+asgn(r ] )r r > (10) 
can be written under the form 

e» = 2 -^(a\v\)^ l^-j , (11) 

where N v depends on the radial n and orbital I quantum numbers, as well as on rj 
a priori. An interesting feature of the AFM is that the functional form of (fTTj) does 
not depend on the chosen potential P{r), provided that this potential is of power-law 
form (typically one can choose r 2 or — 1/r) 3 . N v depends only on the particular 
form of P(r): One logically obtains N v = 2n + 1 + 3/2 or n + l + 1 according to the 
choice P(r) = r 2 or —1/r. Since the form of N v is an artifact of the AFM, it can be 
modified to better fit the exact results. In particular, we have shown that, for any 
physical values of 77 (77 > —2), a good form for is given by 

N v = b(ri)n + l + c(r)). (12) 

Simple expressions for the functions 6(77) and c(r/) have been proposed in [2] in order to 
give an approximation as precise as 10 -3 for the most interesting (the lowest) values 
of the quantum numbers n and I. For example, 

4177 + 86 c()= 5^+17 
K " 1377 + 58' KU 277 + 14' V ' 

leads to very accurate energy formulae [2J. Finally, it is worth mentioning that the 
energy formula (fTTj) together with (|12[) and (fTUj) gives the exact result in the following 
important cases: 

• the harmonic oscillator, for which N2 = 2n + I + 3/2; 

• the Coulomb potential, for which 7V_i —n + l + 1. 

As we have analytical and accurate energy formulae for any power-law potential 
at our disposal, we will assume in this work that P^' (r) = sgn(?7) r v is a good starting 
point to apply the AFM. Setting P(r) — P^^(r), we have shown in [3] that the energy 
spectrum of Hamiltonian 

H=^— + V{r) (14) 



is given by 



2m 



E{v ) = ^-v I{ Vo f + V{I{v )), (15) 
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the optimal value vq being determined from the equation 

N 2 

\ri\v I(v ) ri+2 = (16) 
m 

and the function I(i>) coming from the relation 

HvHyf- 1 = V'{I[y)). (17) 

I(v) depends only on the potential V(r) and not on a particular eigenstate, while 
v>o depends both on the potential (through the I function) and on the state under 
consideration (through N v ). Equations (p~5|) - (P^Tjl are general results of the AFM, that 
will be used intensively in this work. 

We have shown in [3] that replacing the power r\ by A does not change the 
functional form of l{v) but substitutes simply the number N v by N\. The AFM 
cannot say what is the optimal value for the power of P(r). It is necessary to resort 
to a numerical analysis. Moreover it was also shown in [3] that the approximate 
eigenvalues are ruled by the correct scaling laws. 



3. Pure exponential potential 

3.1. Energy spectrum 

The Schrodinger equation with pure exponential potential corresponds to the 
Hamiltonian ((21) in which A = 0, i.e. 

n e =q 2 -ge~ x . (18) 

The negative eigenenergies e(g) of such a Hamiltonian can be analytically computed 
for I = only. But even in this case, the expression of e(g) is not very tractable since 
it is formally defined by the relation 1 

■V=^)(VS) = 0, (19) 

where J p is a Bessel function of the first kind. Consequently, it is of interest to find an 
analytical energy formula for the energy levels of the pure exponential potential that 
is of simpler use than (|19p for I = and, above all, that remains valid for arbitrary n 
and / quantum numbers. 

Equations ([T5)) - (fTT]) can be particularised in the present case. Obviously one 
has to set V(x) = —ge~ x and m — 1/2, while a particularly convenient choice is 
P(x) = PW(ar) = x. A quick calculation leads to 

I( V )=-la(Z) (20) 



and to the energy formula 

-Ur, — , 

2 \g 



e e (*/ ) = -^o-?m(-], (21) 



where vq satisfies the equation 

- 1/ In — 
L V 9 



2N 2 . (22) 



Notice that, in principle, N = iVj in this last equality. But, as mentioned above, 
another choice N\ with an arbitrary value of A is a priori possible. It is worth noting 
that, from ([6]), we have K{x) = ge~ x and thus v € [0,g]. 
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By defining 



x 



1/3 



and Z 



2N< 



1/3 



(23) 



we can recast (|2"2"|) into the form 
— xq In xq = Z. 



(24) 

-xolnxo is such that F(0) 



xq G [0, 1] by definition, and the function F(xq) 
F(l) — 0, with a maximum in F(l/c) = 1/e. It means that (|24|) has no solution if 
Z > 1/e; the solution of (j24|) turns out to be 

x (Z) = c w o(-z) = _ ^ with <Z<l/e, (25) 

where Wo is the principal branch of the Lambert function, whose main properties are 
given in |Appcndix A| The choice of Wq ensures that the total energy has the correct 
limit for small (3 (see below). Equation (|2"Tj) now becomes 



e e (i>o) = e e {gx (Zy 



3W (-Z) 



1 + -W (-Z) 



and the energy, rewritten in physical units, is finally given by 



E„ 



2m 



-ae 



3W (-Z) 



1 + -Wa(-Z) 



with 



N 2 f3 2 



1/3 



(26) 



(27) 



The limit (3 — > is interesting. Indeed, — ae~^ r w —a + a/3r in this case: The 
potential is nearly a linear confining one. Recalling that g — 2ma//3 2 for the pure 
exponential potential, it can be checked that, at small (3, 



E, 



m - 1/3 [a(3N] 



2/3 



0((3 4 / 3 ). 



This last formula is exactly the one which is expected in the case of a linear potential 
[see (fTTj) for a — af3 and A = 1]. 

For computational reasons, the Lambert function may be uneasy to deal with. 
We mention that we can use the approximate expression 

a -3Z 



^3W (-Z) 



1 



\w ( 



-Z) 



100 



[100 -150 Z 



580 Z 2 



524 Z 6 



(29) 



which leads to an absolute error below 0.008 for < Z < 1/e. 

We finally point out that the energy spectrum of potentials like (a — be~ cr ) are 
trivially obtained by adding a constant term to formula (1271) . Such potentials are of 
interest in hadronic physics when one considers a screened confining potential to take 
into account possible string-breaking effects [51 [S] . 



3.2. Critical heights 

Let us consider first the AFM framework. By definition, the critical heights g n i of the 
exponential are such that e e (g = g n i) = 0, with e e given by (|2l)]) . According to this 
last equation, one finds that the energy vanishes for Z = 2c~ 2 / 3 /3, a value which is 
lower than 1/e as required by ([24]) . Equivalently, one can say that the critical heights 
are given by 

9.ynl = \N^. (30) 
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They are such that, if g > g m i , the potential admits a bound state with the given 
quantum numbers n n , Iq. It is remarkable that our approximation scheme, based on 
the linear potential gx for which an infinite number of bound states is present, is 
able to predict that only a finite number of bound states will be present in the pure 
exponential potential. This is another test of the ability of the AFM to analytically 
reproduce the qualitative features of a given eigenvalue problem. Let us note that 
the AFM does not give any information about the optimal power-law potential to 
determine the form of the number N n i . 

Equation (fTi?)) can also lead to a determination of the exact (as coming from the 
Schrodingcr equation) critical heights g n o; let us denote them as <?* . Indeed these 
critical heights are such that they correspond to a energy level e(g^ ) — 0. We are 
thus led to the equation 

Jo(2^~ ) = 0^g* = J -f, (31) 

where j n is the (n + l) th zero of the Bessel function Jo. At large n, these are given by 
j n ps 7r(n + 3/4) 8J, leading to the asymptotic expression 

9* n o *4 (n+l) ■ (32) 



4 



2 

/if)' 



This result is qualitatively similar to (|30p . stating that g n o oc N 1 

We can try to use formula (f32|) to improve the result ((30)) . Assuming that 
N n i = bn + I + c with no constraint on parameters b and c (contrary to what is 
expected from (fT2|) ). it is easy to see that formulae ([30]) and (|32f coincide for g n a with 
b = 7r/e and c = 37r/(4e). We can then try a new relation to compute the critical 
heights 

/ 7r e 37r\ 2 

ge;nl = I ^ n + 3 + J ' ( 33 ^ 

This formula is in good agreement with exact results: For n S [0, 5] and I £ [0, 5], the 
minimal, maximal and mean relative errors are respectively 0.03%, 8.7% and 3.5%. It 
is not possible to obtain so good agreement by choosing b and c values coming from 
(|13| . Some exact g* e . n i values are given in table [TJ they are computed numerically by 
solving the Schrodinger equation with a vanishing energy. 



Table 1. Some exact critical heights g*. nI for the exponential potential — ge x . 



n 


I = 


1 


2 


3 





1.446 


7.049 


16.313 


29.258 


1 


7.618 


16.921 


29.880 


46.518 


2 


18.722 


31.526 


48.077 


68.346 


3 


34.760 


50.948 


71.002 


94.838 



Formula ([33]) is fully analytical and is designed in order to reproduce correctly the 
behaviour of the exact critical heights for large n and / = 0, g e -,n0 ~ ffno- Unfortunately, 
the analytical asymptotic expression of g^ t is not known, and an improvement of 
formula (|33() must rely on numerical arguments, in particular using the numerical 
values of g*. n i, as given in table [T] Moreover, a good overall fit formula needs to 
reproduce correctly not only the behaviour in terms of n or / separately, but must 
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try to get the good behaviour in terms of n and I simultaneously. Since we want to 
stick as much as possible to the philosophy of AFM, we maintain the quadratic form 
g e - n i cx ./V 2 and just add in the expression of N a term proportional to y/rd. Thus, we 
propose a new expression for the critical heights which reads 

g e . nl w (1.566 n + 1.393 1 - 0.125V^7 + 1.202) 2 . (34) 

The improvement as compared to (|33[) is very appreciable since, in this case, the 
minimal, maximal and mean relative errors are respectively 0.005%, 4.5% and 1.1%. 
Using (j30|) , one can write 

Z = |c- 2 / 3 (^) . (35) 

This makes clear the dependence on the critical heights g e -ni for formula (f26|) . It is 
then easy to see that e e = when g = g e ;ni- 

4. The general case 

In order to compute analytically the spectrum of Hamiltonian |J3J) for arbitrary values 
of A, it is of interest to use the general result (HJ), and consequently to choose 
P(x) — sgn(A) x x . Equation ^ can be rewritten in the present case as 

v = K{x)=g^-{x-\). (36) 

The function K{x) is such that K(0) — — gsgn(A), K(X) — 0, and lim^oo K(x) = + . 
It has a maximum in K(X + 1) = </e -1 /|A|, the position of this maximum being in 
the range x > for A > — 1. We can conclude from this discussion that K{x) £ [0, <?] 
if -2 < A < -1, K(x) £ [0,.ge- A -V|A|] if -1 < A < and K(x) e [-g, ge'^ 1 /\X\] 
if A > 0. In the first case, the function K(x) is monotonically decreasing while in the 
two last cases, it presents a maximum. 

The inversion of (|36[) can be achieved thanks to the Lambert function since K (x) 
is of the form (|A.3|) . One finds 



/ A M = A-T4^-c A |A|-J , (37) 

where one must be careful in the choice of either Wo or W-i, that will depend on the 
particular value of A. Nevertheless, as both branches of the Lambert function share 
the same properties, it is possible to go a step further in our computations by using 
the generic solution (|3T|) . 
Let us first define 

|Ahe A 2iV 2 e A 
xq = and Y\ = , (38) 

9 9 
where N — N\ a priori. But let us recall that another choice is possible as mentioned 
previously. Then, a rewriting of (p~5|) - (jTT)) with V(x) = —gx x e~ x , m = 1/2, and 
P(x) = P w (x) = -sgn(x)x A leads to 
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where Xq is a solution of 

x [X-W(-x )} X+2 =Y X . (40) 
When the function Fx(x ) = [•*• — W {— x )] x+2 has a maximum, it is located 



at 



1 



x x = a x e~ a \ with a x = -(V9 + 4A + 3). (41) 



Notice that (3 — V9 + 4A)/2 would also be a possible value for a\, but the choice (|4"Tj) 
ensures that 2 a is located in the domain of Xq. One finds the constraint 

Y X <F X = F{x x ) = a x e- a > (A + a x ) X+2 . (42) 

As for the exponential potential, (l39|) vanishes for particular values of Y x , and 
thus critical heights appear. It can be checked that ex is equal to zero when 
Y x = 2(A + 2) A+2 /e 2 , which is lower than Fx as demanded by (|4"2"]l. The critical 
heights are then given by 

/ e \ A+2 

9x-,m = [xT2j N " 1 - (43) 

Again, the AFM gives no information about the best form for the number N n i. A first 
trial is to use N — Nx but other choices are possible. The case A = was discussed 
in the previous section. The case A = — 1 will be studied in the next section. 

As a general result, (|43|) states that g n i oc for any exponential potential. 
Interestingly, it has previously been observed that the critical heights of an arbitrary 
exponential potential can be very accurately fitted by a biquadratic form in n and / [9] . 
Our results can be seen as a demonstration of this empirical observation. However, 
the accuracy of the fitted form in [9] is a supplementary argument showing that the 
explicit form of N n i cannot be computed within the AFM: Only N n i = bn + I + c can 
be assumed, the coefficients b and c having to be determined by another procedure in 
order to improve the analytical formulae that we have obtained. 

Equations (f3"9"]l and (|40p formally define the analytical energy formula we are 
looking for. Finding a general solution of (|40|) is a complicated problem that cannot 
be solved analytically. Moreover, the Yukawa potential is, to our knowledge, the only 
potential of common use with A ^ 0. Before making a full treatment of the case 
A = — 1 in the next section, we can nevertheless make some comment about the limit 
where /3 goes to zero (g — > oo and Y\ — * 0). In this limit, —a r A e~^ r as —a r x , and we 
should recover Ijlip for A < and a — a. The nontrivial solution of ([40]) for Yx = 
and A < is Xo = — Ae A . It can then be checked that 

o2 9 , \ / N 2\ VO+ 2 ) 

km E x = Mm ^£a(-Ac a ) - l±^( a |A|) 2 /( A + 2 ) ( — ) .(44) 

This is the expected result for the energy spectrum of the potential — a/r' A l. It is 
important to notice that the trivial solution xq = would lead to the unphysical limit 
Ex -» 0. 

In the limit A —> 0, equation (|40|) reduces to 

x W(~x ) 2 =Y Q = (3Z) 3 (45) 
and its solution is given by 

Xo = ~3W(-Z)e 3W{ - z \ (46) 
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in which Z is defined by 
Taking into account that - 
reduces to 



(paying no attention how the number N is defined). 
Z e [-l/e,0], a range in which W(-Z) < 0, and that ([55]) 



(47) 



one recovers formula (|26p . but with a priori N — N . This is an illustration that 
the same functional form is obtained independently of the power-law chosen for the 
potential P(x). Moreover, g\ = o-jii = g e -ni and (|5D|) is also recovered. 

Up to now, the branch to be used in these formulae is not yet determined. For 
— 2 < A < — 1, the function K (x) defined by ([36]) is monotonic and its inverse is single 
valued; a careful inspection of its properties shows that only the W-\ branch must 
be retained. For A > — 1, the situation is more complex and we did not find a simple 
prescription to fix unambiguously the type of branch. Fortunately the corresponding 
potentials are not of common use. 



5. Yukawa potential 

The Yukawa potential corresponds to the case A = — 1. A look at ([36]) in this case 
shows that one has vq € [0, g] and equivalently xq €E [0, 1/e]. Since x = I-i(v) £ [0, oo[, 
W-\ must be chosen: This branch of the Lambert function is defined on a finite 
domain, but can go to infinity as demanded. In this case, eqs. ([39]) and ([40]) read 

ey{Y) - ~ g 2Y[xo(Y) + Y] ' (48) 

- x Q (Y) [1 + W-x(-xo{Y))] = Y, (49) 

where Y denotes Y-i = 2N 2 /(eg). Notice that, in principle, N = N-i in this last 
equality. But, as mentioned above, another choice N n with an arbitrary value of 
77 is a priori possible. The function F-i(x ) = —x [1 + W-i(— x )] is such that 
.F_i(0) = .F_i(l/e) = 0, with a maximum in 

^ = ^e-¥, F = F_ 1 (x) = (2 + V5)e-^. (50) 

Let us note the relation F/x = $ = (1 + \/~5)/2, the golden ratio. Following (|43|) . the 
critical heights read in this case 

g v -ni = eN^. (51) 

A previous formula giving analytical approximate energy levels for the Yukawa 
potential has been found in [lOj in the framework of the envelope theory. It is given 
in our notations by 



ll x>a \ x z " x 



9 • (52) 



But, we show in |Appcndix B| that E^f is exactly equal to e y (Y) defined by (05]) and 
(|49[) . It means that the AFM and the envelope theory lead to the same approximate 
energy formula for the Yukawa potential, provided the quantity N is restricted to 
N-\. Nevertheless, a closed energy formula has not been proposed in [TO]; we achieve 
such a task in the following. 
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Figure 1. Plot of the exact solution of (149 I t (solid line), compared to the fitted 
form l(53)l with ^4 = (full gray line), A = 2 (dotted gray line), and A = A c ss 2.87 
(dashed gray line). The case A = 2 is nearly indistinguishable from the exact 
solution. 



In order to have an explicit energy formula, it is of interest to find an analytical 
approximation for xq(Y), the solution of (|49f . For a given value of Y, two values of Xq 
are allowed: The first one goes to zero when Y vanishes, and the second one tends to 
1/e in this limit. We already mentioned that the physical solution should actually be 
such that lim-p^Q Xq(Y) — 1/e in order to find the correct Coulomb limit when (3 — > 0. 
From the study of the function F-i(x), it is possible to determine that xo(0) = 1/e 
and xq(F) = x. Moreover, in Y = and Y = F, the function xq(Y) has a horizontal 
and a vertical tangent respectively. So, the simplest form to parameterise xq(Y) is a 
branch of ellipse. 

As it can be seen in figure [TJ the analytical form 



x (Y,A) = x + 




AY(Y — F) 



(53) 



reproduces rather well the numerical solution for the simplest choice A = 0, but is 
particularly accurate for the optimal value A = 2. When 4^0, the curve is not a 
branch of ellipse and the tangent in Y = is not horizontal, but the global agreement 
is improved. Recalling that g = 2ma/ (3 for the Yukawa potential, the final energy 
spectrum reads 



E y = ^e y (x (Y,2)) 



-a(3 



ex (Y,2) 2 [x (Y,2)-Y} 



2Y[x Q (Y,2)+Y] 
This formula can be recast in the following form 

. 2 x (Y,2)-Y 



with Y = 



I3N 2 



e ma 



E„ — E, 



Coul 



e 2 x (Y,2Y 



x a (Y,2) + Y_ 

where -Ecoui is the energy for a pure Coulomb potential given by 



a 2 m 
'2N 2 ' 



(54) 



(55) 



(56) 
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As the term within brackets in {55J tends toward 1 when (3 — > 0, we can expect that 
N = 7V_i is a good choice for small values of (3. This will be examined in section [5] 
Using (|5ip , one can write 

y = 4— ■ (57) 

e 2 3 

This makes clear the dependence on the critical heights g y ^ n i for formulae (j5~4")l and (f55|) . 
It is worth mentioning that A = 2 comes from a fit of formula (|53[) on the numerical 
solution of (|49p . Another choice could have been done, that is to set A = A c such that 

z Q (y = ^,A C ) = %. (58) 



. 2 ' V e 2 ' 

Then the final energy formula exactly vanishes at the critical heights (fBTj) . One then 
obtains for A c a complicated analytical formula that we do not give explicitly here; it 
is enough for our purpose to know that A c ss 2.87. 

The results we obtained here in the most interesting case of the Yukawa potential 
can actually be generalised to potentials with — 2 < A < — 1. It can indeed be 
numerically checked that the solution of (|40[) in this range is quite well fitted by 
the analytical form 



x (Y x , A x ) = x x + (|A|e A -x x )Jl- (j^ 



A X Y X (Y X -F X ), (59) 



with 

A x = —(109 + 196A + 85A 2 ). (60) 

This function is such that A-± — 2. Notice that the branch W-i has to be chosen, 
by continuity with the Yukawa case in particular. Once this solution is assumed, the 
energy spectrum is eventually given by e x (xo(Y x , A x j) [see (|39| with W = W-i]. 

It has been shown in [11 that the critical heights of the Yukawa potential are 
given with a relative accuracy around 0.4% by the following fitted form 

\ 2 

f 

where 



«^=2(v^ + -) , (61) 



Zi = Z (l + al + f3l 2 ), Si = S (l + -fl + SI 2 ), 
Z = 0.839908, a = 2.7359, 0= 1.6242, 

So = 1.1335, 7 = 0.019102, 6 = -0.001684. (62) 

Due to the smallness of 7 and 5, we can replace Si by Sq in (fBTjl . It clearly appears 
that g n i oc (iSoV Zq/3 I + n) 2 asymptotically. This quadratic behaviour is predicted by 
our analytical results. Introducing a constant which allows to obtain the exact result 
for n = I = and considering the asymptotic expansion for large values of I, (|61[) is 
approximately given by 

g% w (1.248 n+ 1.652 Z + 1.296) 2 . (63) 

It could be interesting to improve the values of critical heights for the Yukawa 
potential as it was made for the pure exponential potential. Unfortunately, no 
analytical formula for g y - n o is available in this case. Nevertheless, we can assume 
that the formula (fSTj) is a good starting point and that N n i = bn + I + c with no 
constraint on parameters b and c (contrary to what is expected from (jX2j) ) . If the 
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two exact values gy. 00 and g^-io are known, the choice b — y5y-io/ e — y9y-oo/ e an< ^ 

c = yVy ; oo/ e msures that (jSTj) gives the correct critical heights for I = and n = 
and 1. Instead of using exact numerically computed values, we can use analytical 
approximate evaluations of these quantities (see | Appendix CJ ). We can then try a new 
relation to compute the critical heights 

9y,m = ([V^vixo - y/gy,oo] n + y/el + \/^oo) ■ (64) 
With the exact critical heights, we have 

g y . n i « (1.243 71+ 1.649 Z+ 1.296) 2 , (65) 
while with the use of approximations (|C.2[) . we have 

g y;nl w (1.291 n + 1.649 1 + 1.296) 2 . (66) 

These formulae are in good agreement with (f6"3")l and with exact results: For n € [0, 4] 
and Z 6 [0, 4], the minimal, maximal and mean relative errors are respectively 0.005%, 
4.5% and 1.2% with §Q and respectively 0.04%, 5.8% and 3.4% with §6$. It is not 
possible to obtain so good agreement by choosing b and c values coming from (I13[) . Let 
us note that the minimal, maximal and mean relative errors are respectively 0.01%, 
0.9% and 0.3% with formula ([61]) . The very good quality of these results is due to the 
use of a complicated formula fitted on exact results. Our formulae are simpler and 
the general behavior is predicted by the AFM. Some exact g*. nl values are given in 
table [U they are computed numerically by solving the Schrodinger equation with a 
vanishing energy 



Table 2. Some exact critical heights g*. nl for the Yukawa potential —ge x /x. 



n 


I = 


1 


2 


3 





1.680 


9.082 


21.895 


40.136 


1 


6.447 


17.745 


34.420 


56.514 


2 


14.342 


29.461 


49.970 


75.899 


3 


25.372 


44.261 


68.572 


98.318 



In the very same way that we improved the critical heights of the pure exponential 
potential, it is possible to improve the corresponding Yukawa quantities with the same 
kind of prescription. Here we propose the expression 

g y;n i~ (l.247n+ 1.680Z-0.054V^7+ 1.296) 2 . (67) 

In this case, the minimal, maximal and mean relative errors are respectively 0.002%, 
2.4% and 0.6%. We obtain a quality of the same order than the one obtained with 
formula (|61[) . but with a much simpler expression. 

6. Improved energy formulae 

In the previous sections, approximate analytical forms for solutions of a pure 
exponential or a Yukawa potentials were found. The formulae depend on the quantum 
numbers n and I through a factor N. This number could be taken as N\ without 
information about the optimal value for A. If we look for instance at the Yukawa 
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potential, it is clear that it reduces to a pure Coulomb interaction when /3 = 0. In 
this case, the choice N = iV_i gives the exact result in this limit. 

The dimensionless Hamiltonians considered above depend on a parameter g. The 
variation of the eigenvalues being smooth for the variation of g, we can assume that 
the number N, giving the optimal values for all the eigenvalues of a Hamiltonian, is 
also a smooth function of g. From considerations above, the functional form 

N(g) = b(g)n + l + c(g) (68) 

seems reasonable. As we cannot predict the correct behaviour, it is necessary to focus 
our attention on numerical solutions. 

Very accurate eigenvalues enum(ff; n i ^) for Hamiltonians defined above can be 
obtained numerically with the Lagrange mesh method [12j . It is very accurate and 
easy to implement. In order to find the best possible values for coefficients d{g) (d 
stands for b or c), we will use the following measure 

X(9) = ^2 ( £ num(.g; n, t) - eapp(g; n, t)) , (69) 

{n,l} 

where £app(ff! n,i) are values obtained from our approximate formulae. Other choices 
are possible but we find this one very convenient. The summation in (|69[) runs a priori 
on all possible pairs {n, 1} for the bound states allowed by the value of g. But for some 
special pairs {ro, I}, the value £app(<?; ti, £) can be a complex number (condition (|40p is 
not fulfilled) or a positive number (quality of approximation is poor). This drawback 
can happen for the highest levels. In this case, the pair {n, /} is not taken into account 
in @g). 

The analytical form eapp depends on N(g) which depends on coefficients d. For 
each value of g, optimal values for the d coefficients, d m ^ n (g), can be determined by 
minimizing x(<?). Then, with a set {^nrm(<?)} for a given set {<?}, a functional form 
dfu(g) can be fitted with the following measure (%' is also a chi-square fit, but the 
prime indicates that the function \' is indeed different from \ i n (|69|) ) 

x'(d) = J2(d min (g)-d flt (g)) 2 . (70) 

{g} 

Again, other choices are possible but we find this one very convenient. We will now 
try to determine the best form of coefficients d(g) for the potentials studied above. 



6.1. Pure exponential potential 

Approximate dimensionless eigenvalues for the pure exponential potential are given 
by ([26)) . By minimising our measure x{9)i we found the optimal values of b{g) and 
c{g) for several finite values of g. The results are plotted with dots on figure O The 
small irregularities in the data are due to the fact that the number of bound states is 
obviously not a smooth function of g. 

We tried to fit the numerical points with various functions and found that 
the best result is obtained for sections of hyperbola, for both b(g) and c(g). The 
parameterisation retained is 

1.42g- 12.76 . . 1.32 g+ 16.88 
b{9)= ff -8.62 ' C{9)= 5 + 14.95 ' 

The agreement with numerical points is good as it can be seen on figure [2j For large 
value of g, a constant value of the coefficient b is a good approximation. This is clearly 
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not the case for the coefficient c. The quality of this parameterisation for both 6(3) 
and c(g) in can be appraised in tables [3] and SJ Not all bound states can be found. 
As in (|69p , only real strictly negative eigenvalues are retained. The mean relative error 
is about 20% but the accuracy is much better for the lowest eigenvalues. 




20 40 60 80 100 20 40 60 80 



9 9 

Figure 2. Best values of the coefficients b(g) and c(g) to parameterise the 
eigenvalues H26H for a pure exponential potential: numerical fit with II69I I (dots); 
functions l|71[ l (solid line). 



Table 3. Quality of various approximations for a pure exponential potential. 
R: (number of negative eigenvalues found)/ (total number of bound states); 
"^min /max : mm i ma V maxrma l relative error (%) for negative eigenvalues found; 
Amean: mean relative error (%) computed with all negative eigenvalues found. 
First line: formula l|26| l with parameterisation <|71fl ; Second line: formula 1(26} 
with parameterisation 11721 . 



9 


5 


10 


20 


40 


60 


80 


100 


R 


1/1 


3/3 


4/6 


7/10 


12/15 


15/20 


19/26 




1/1 


3/3 


4/6 


8/10 


12/15 


16/20 


21/26 


^min. 


0.8 


4.1 


1.7 


1.6 


0.4 


0.5 


0.3 




2.0 


1.8 


0.05 


0.3 


0.07 


0.09 


0.2 


A max. 


0.8 


54.6 


79.7 


35.0 


98.5 


91.2 


81.7 




2.0 


33.4 


6.2 


31.7 


80.1 


69.1 


97.4 


Amean 


0.8 


31.6 


22.6 


11.8 


23.6 


18.3 


19.1 




2.0 


14.0 


2.4 


8.9 


11.8 


11.2 


17.2 



By using the form (|68p . we restrict the number of functions to fit and we keep the 
same level of numerical complexity as the one used in [3]. Since relation (|34[) works 
very well to predict the critical heights, we used the same form to fit the eigenvalues. 
Performing again similar calculations as above with (|69|) and (|70j) . we found that a 
good parameterisation is given by 

. . 1.445- 11.17 0.95g- 1.36, 1.43.g + 23.09 
N{s) = .9-6.86 " + 3-0.33 1 + 3 + 20.60 

+ (6.69 10~V - 0.00019 3 - 0.126)\/n7. (72) 

Its quality can be appraised in tables [3] and |4] The eigenvalues are generally improved 
and the number of bound states found is higher. Note that the influence of the term 
\fnl is small. The neglect of this term does not spoil very much the results. 

When — > 0, that is 3 — > 00, the pure exponential potential reduces at first order 
to a constant and at the next order to a linear potential. Relations (f7Tj) show that 
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linig^oo b(g) = 1.42 and that linig^oo c{g) = 1.32. These numbers are different from 
the values b w 1.79 and c w 1.38 expected for a linear potential, and from b w 1.48 
and c w 1.21 expected for constant potential (see l|13p). We will see that the situation 
will be clearer for the Yukawa potential. 

Table 4. Eigenvalues for a pure exponential potential with g = 40 as a 
function of {I, n} sets. First line: exact value; Second line: formula 11261 with 
parameterisation 171t ; Third line: formula (1261 with parameterisation 11721 . A * 
indicates a non real or a non negative value. 



/ 


n = 


1 


2 


3 





-17.53 


-6.88 


-1.87 


-0.077 




-17.84 


-7.58 


-1.59 


* 




-17.30 


-7.28 


-1.43 


* 


1 


-10.14 


-3.35 


-0.42 






-9.98 


-2.99 


* 






-10.10 


-3.55 


* 




2 


-5.03 


-0.93 








-4.63 


* 








-5.06 


-0.64 






3 


-1.55 










-1.01 










-1.52 









6.2. Yukawa potential 

Approximate dimensionless eigenvalues for the Yukawa potential are given by (|48|) 
supplemented by (f53|) . In all results presented, we chose A = 2. By using a similar 
procedure as those described for the pure exponential potential, we found that the 
best result in fitting the numerical points is obtained for sections of hyperbola, for 
both b(g) and c(g). The parameterisation retained is 

0.99 g- 5.92 1.00 g- 1.68 

g — 5.08 g — 1.58 

The agreement with numerical points is good as it can be seen on figure[3] As expected, 
N(g) — > N-i when g — > oo, that is to say when (3 — > 0, since in this limit the 
Yukawa potential tends toward a pure Coulomb interaction. So, in the following, we 
will also consider the choice N = N-i, that is b = c = 1. Again a more complex 
parameterisation 

, r/ , 0.99o- 7.16 1.00 + 2.51, 1.00 0-1.89 

N(q) — n H 1 -\ 

yy ' + 3.16 .9- 1.79 

+ (-0.000233 g 2 + 0.0202 g ~ 0.480) V^T. (74) 

is also studied. 

An accurate energy formula has been found in [TT] from a fit of the numerically 
computed energy levels of the Yukawa potential. In our notations, it reads 

where A' = 1.9875, B' = 1.2464 and a = 0.003951, where N ni = n + I + 1 as in the 
Coulomb case, and where o5 is defined by (|6ip . This formula is rather different from 
ours, given by (|48|) . However, they coincide at the limit g — > oo. 
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9 9 

Figure 3. Best values of the coefficients b(g) and c(g) to parameterise the 
eigenvalues l|48| l with 153|l for a Yukawa potential: numerical fit with 169|l (dots); 
functions (1 73 I t (solid line). 



Table 5. Quality of various approximations for Yukawa potential. R: (number 
of negative eigenvalues found)/ (total number of bound states); A m j n y max : 

minimal/maximal relative error (%) for negative eigenvalues found; Amean: mean 
relative error (%) computed with all negative eigenvalues found. First line: 
formula 1481 and 1531 1 with parameterisation 1731 ; Second line: formula 1481 1 and 
153( with N = N—i; Third line: formula l|48J l and l|53| l with parameterisation (|74j l: 
Fourth line: formula 1751 1. The abnormally large values marked by t are due to 
the existence of a nearly zero eigenvalue which is badly reproduced. 
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5 


10 


20 


30 


40 


50 


R 


1/1 


2/3 


3/5 


6/8 


6/10 


10/13 




1/1 


1/3 


3/5 


6/8 


6/10 


10/13 




1/1 


2/3 


3/5 


6/8 


6/10 


10/13 




1/1 


3/3 


5/5 


8/8 


10/10 


12/13 


^min. 


0.2 


0.07 


0.0006 


0.02 


0.02 


0.02 




13.5 


3.5 


1.3 


0.8 


0.6 


0.5 




0.7 


0.2 


0.01 


0.01 


0.02 


0.02 




0.4 


0.1 


0.06 


0.04 


0.03 


0.02 


A max. 


0.2 


4.3 


6.6 


27.8 


9.4 


42.4 




13.5 


3.5 


14.1 


45.6 


19.2 


57.3 




0.7 


8.6 


0.7 


30.0 


9.8 


44.1 




0.4 


10.9 


14.0 


34.9 


575. 2t 


24.6 


Amean 


0.2 


2.2 


3.2 


13.0 


4.6 


16.2 




13.5 


3.5 


7.7 


19.6 


8.2 


21.2 




0.7 


4.4 


0.3 


6.8 


2.1 


11.1 




0.4 


6.4 


3.3 


9.8 


59.5 1 " 


3.9 



The quality of the various approximations can be appraised in tables [5] and [5] 
Not all bound states can be found, even if formula (|75|) gives generally the correct 
number of bound states. As in (|69j). only real strictly negative eigenvalues are retained. 
Formula (|48|) with approximation (|53|) and parameterisation ([73]) gives quite good 
results. The mean relative error is about 10% but the accuracy is much better for 
the lowest eigenvalues. The parameterisation (|74p gives generally better results. The 
quality of the fit for formula (|75[) and our approach is in average quite comparable. 
But formula (f75|) is empirical. Conversely, (|48|) is obtained from an explicit analytical 
resolution of the Schrodinger equation. Our results are not so good by choosing the 
parameterisation TV = iV_i but the quality is still reasonable. 
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Table 6. Eigenvalues for a Yukawa potential with g = 30 as a function of 
{I, n} sets. First line: exact value. Second line: formula 148H and d 53 1 ) with 
parameterisation J73}. Third line: formula (1481 and 153H with N = iV_i; Fourth 
line: formula 1481 and 1531 1 with parameterisation 174H ; Fifth line: formula 1751 1. 
A * indicates a non real or a non negative value. 



i 


n = 


1 


2 


3 





-196.44 


-31.51 


-5.47 


-0.22 




-196.41 


-32.31 


-4.16 


* 




-194.82 


-29.66 


-2.97 


* 




-196.42 


-31.74 


-3.89 


* 




-196.36 


-31.62 


-5.67 


-0.30 


1 


-30.74 


-4.94 


-0.029 






-29.84 


-3.57 


* 






-29.66 


-2.97 


* 






-30.89 


-4.82 


* 






-30.52 


-4.85 


-0.019 




2 


-3.81 










-3.02 










-2.97 










-3.49 










-3.70 









7. Summary of the results 

The AFM was proposed in [2 , 3J as a tool to compute approximate analytical solutions 
of the Schrodinger equation. The basic idea underlying this method is to replace an 
arbitrary potential V(r), for which no analytical spectrum is known, by an expression 
of the type v P(r) +g(v) 1 P(r) being a potential for which analytical eigenenergies can 
be found, v the auxiliary field and g{y) a well-defined function of this extra parameter. 
This auxiliary field is such that its elimination as an operator leads to the original 
Hamiltonian. If v is seen as a number however, analytical approximate solutions can 
be found in favourable cases, and the auxiliary field is eventually eliminated by a 
minimisation on the eigenenergies. It is shown in [21 [3] that this method is a kind of 
mean field approximation. 

In the present work, we have searched for solutions of the Schrodinger equation 
with exponential potential of the form —ar x e~ l3r using the AFM. General analytical 
results have been obtained for the eigenvalues of this type of potentials. It is 
proved that the critical heights for these interactions are proportional, with a good 
approximation, to (bn + I + c) 2 where b and c are constant parameters. 

We have focused our attention on pure exponential and Yukawa interactions. 
Approximate analytical formulae have been found for these two particular cases for 
arbitrary values of n and I. Some results were previously obtained for the Yukawa 
potential [101 [13] but, up to our knowledge, it is the first time that a closed energy 
formula is computed from the Schrodinger equation. The majority of eigenvalues can 
be computed. The accuracy is quite good, mainly for lowest states. Moreover, very 
simple analytical approximate formulae have been obtained for the critical heights of 
these two potentials. 

As an outlook, we mention that more complicated problems could be studied with 
the AFM. In particular, we plan to apply this method to semirelativistic Hamiltonians. 
Such a work is in progress. Moreover, we have shown that the energy spectrum we get 



Auxiliary field method for exponential potentials 



19 



with the AFM for the Yukawa potential is exactly equal to the one which is obtained 
with the envelope theory. This suggests the existence of a connection between both 
methods, that we will discuss in a subsequent paper. 
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Appendix A. The Lambert function 

Let us briefly recall some points concerning the Lambert function (also called Omega 
function or product-log), that we will denote W(z). First of all, W(z) is defined as 
the inverse function of ze z . Consequently, it has the following properties: 

W{z)e w(z) = W{ze z ) = z, (A.l) 

8 - m! » = iiirli' < A - 2 > 

But, it is readily observed that the inverse function of z e z is multivalued. Two 
branches of the Lambert function thus exist, respectively denoted as Wo(z), defined 
for z > — 1/e, and W-i(z), defined for — 1/e < z < o[J They are plotted in figure [KT\ 
Obviously, both branches of the Lambert function share the same properties (|A.1[) 
and (|A.2|) as W(z); they meet in W-i{— 1/e) = Wo(—l/e) = — 1. It can moreover be 
checked that Wo(|x| <C 1) ~ x, lim^o W-i{x) — — oo, and lim^oo Wq(x) = oo. 



X Notice that the two branches Wq(x) and W—i(x) of the Lambert function W(x) are known by the 
software Mathematica package as ProductLog[0,x] and ProductLog[-l ,x] respectively. 
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Figure Al. Plot of the two branches of the Lambert function, namely Wo(z) 
(solid line) and W—i(z) (dotted line). 



For our purpose, it is worth mentioning that the equation 
(az + b) n c- z = 
is analytically solvable with the Lambert function. One has 



z= nW 

a 



.JL( e -V«0)iA. 
a n 



(A.3) 



(A.4) 



where either Wo or W-\ has to be chosen following the range of 9 and z that the 
function z{9) has to cover. 



Appendix B. Yukawa potential and envelope theory 

The envelope theory is a powerful method allowing to get upper and lower bounds on 
the energy levels of a given potential [H] . The case of the Yukawa potential has been 
studied in [10]. In particular, it is shown that an upper bound of the energy levels is 
given in our notations by 



E 



up _ 



f(x) with f(x) = — g . 

x 



an>0 ' • x 

Equivalcntly, one can write the energy formula as 

Kf = f( & ) with /'(0) = o. 

Explicitly, these relations can be rewritten as 

nl 6 2 l 1 + 



1 e" 



Y. 



(B.l) 
(B.2) 

(B.3) 
(B.4) 
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These last two equations are completely equivalent to (H51) and P5|) , giving the 
energy spectrum of the Yukawa potential obtained with the AFM, provided that the 
following redefinition is made 

x = (e + l)c- ( - 9+1 \ (B.5) 

and that the quantity iV_i is replaced by N. 



Appendix C. Analytical estimation of critical heights g„o f° r the Yukawa 
potential 

In [15j . it is shown that, for a fixed value of the parameter a, solutions of the equation 

d 2 (/> ( , e" 
dx 2 



^+(* + b— )* = 0, (C.l) 



with the boundary conditions 0(0) ~ 0(oo) = 0, can be accurately obtained by a 
variational method. The eigenvalues b are given by the extremum of the ratio J/N, 
where the quantities J and N are respectively given by equations (18) and (19) in [15] 
(note a misprint: parameter b must not be present in (19)). This ratio J/N depends 
on a and onn+1 linear variational parameters h v . The critical heights g n o we search 
for are given by lim Q ^o b. 

For n = 0, the only solution, which is independent of Iiq (normalisation condition), 
is given by goo = l/ln(16/9) ~ 1.738. This value if not too far from the exact one 
(see tabled]). For n = 1, two extrema for J/N exist and depend on hi/h . They can 
be obtained analytically, but their expressions are very complicated and unusable in 
practice. Nevertheless, the two optimal values for hi/ho are very close to —2 and 1/2 
respectively. So we can set 

J 

N 



9oo 



17 

1.680, 



.9io 



J 

N 



fcl /ho=V2 61n(27/5) 

' 6.693. (C.2) 



hl /h =-2 601n(2)-261n(3)-81n(5) 
These two values are good approximations of the exact ones. 



